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TECHNICAL NOTE 4008 


A SUMMARY OF GROUND-LOADS STATISTICS 

By John R. Westfall, Benjamin Mllwitzky, 
Norman S. Silsby, and Robert C. Dreher 


SUMMARY 


This paper briefly summarizes the more important statistical data 
obtained by the NACA on the subject of ground loads. The information 
presented relates primarily to landing-impact and taxiing loads; however, 
some limited data are also presented on one phase of gr mmfl - hanrn i ng 
loads, namely, braking friction. A number of experimental and theoreti- 
cal papers dealing with various aspects of the subject are listed in the 
bibliography. 


INTRODUCTION 


The NACA has for seme time been concerned with the study of ground 
loads on aircraft; these include landing-impact loads, taxiing loads, and 
loads due to ground handling. A number of experimental and theoretical 
reports that have been published on various aspects of the subject are 
listed in the bibliography. The present paper summarizes the more impor- 
tant statistical data on ground loads obtained by the NACA, which may 
serve as a basis for predicting the ground-loads experience of aircraft. 

The two major phases of the ground-loads problem considered herein 
are landing impact and taxiing loads; however, some limit ed data are also 
presented on one phase of ground-handling loads, namely, braking friction. 

The section on landing impact considers first the initial contact 
conditions, then landing-gear reactions and airplane response. The con- 
tact conditions are in the form of statistical data on five par ame ters ; 
vertical velocity, horizontal velocity, bank angle, rolling velocity, 
and wing lift. Although there are other parameters, the ones considered 
appear to be the most important. The discussion of ground reactions anrj 
airplane response is primarily concerned with the analytical prediction 
of vertical and drag loads from the known initial conditions. 

On the subject of taxiing loads, data are available in the form of 
runway profiles and acceleration measurements from VGH records in taxi ing ; 
the latter are presented in statistical form. Some considerations on air- 
plane structural response in taxiing are also given. 
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With regard to ground handling, some limited data are presented on 
the coefficient of friction as measured in braking tests with an airplane 
and with a specially designed tow cart. 
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SYMBOLS 

vertical ground force, lb 
frequency, cps 

acceleration due to gravity, ft /sec 2 

incremental center-of -gravity acceleration, g units 

tire pressure, lb/sq in. 

vertical velocity at contact , fps 
airplane weight, lb 
coefficient of friction 


Subscript : 

max maximum 


RESULTS AND DISCUSSION 
Landing Impact 


Statistics.- The available statistics on the initial, contact condi- 
tions are summarized in figures 1 to 4. Figure 1 shows probability curves 
for the vertical velocity at contact. The ordinate is the probability of 
equaling or exceeding a given value of vertical velocity; the abscissa is 
the vertical velocity. Because there appeared to be definite differences 
in the statistics for civil and military (land-based) airplanes, the prob- 
ability curves for these two categories are shown separately. Although 
the reasons for the differences between the two categories are not com- 
pletely understood, it is probably significant that, on the whole, two 
different generations of airplanes are involved - the civil airplanes 
were all piston powered, whereas the military airplanes were mostly jets. 
Furthermore, some of the military landings were transitional-training 
flights. The more severe character of the military landings is indicated 



NACA TN 4008 


3 


by the generally higher values of vertical velocity at a given probability 
level. For example, 1 in 10,000 landings of military airplanes would be 

expected to equal or exceed a vertical velocity of about 8^ feet per sec- 
ond; the corresponding value for the civil airplanes is about 5 feet per 
second . 

Figure 2 shows s imila r probability curves for the airspeed at con- 
tact. The ordinate is the probability of equaling or exceeding a given 
value of airspeed; the abscissa is the airspeed at contact, expressed as 
percent above the stalling speed. Again, the military operations appear 
to give higher velocities at a given probability level. 

With regard to bank angle and rolling velocity at contact, which are 
of concern for unsymmetrical landings, wing span appeared to be a better 
criterion for separating the data than types of operations. Accordingly, 
with the number of engines taken as a crude measure of wing span, the 
airplanes with four and more engines were put into one category and the 
single- and twin-engine airplanes into another. Figure 3 presents the 
probability curves for bank angle. As might be expected from geometric 
considerations, the larger span airplanes have somewhat smaller angles 
of bank. The probability curves for roll ling velocity are shown in fig- 
ure 4; here the differences are more pronounced than for bank angle. It 
may be that these differences reflect not only the geometrical span 
effect, but also the greater moments of Inertia and damping in roll of 
the larger airplanes. 

A study was made to determine whether there was any statistical 
correlation between bank angle and rolling velocity, on the one hand, 
and vertical velocity or horizontal velocity, on the other hand. It was 
found that no such correlation existed. There was also no correlation 
between bank angle and rolling velocity. It therefore appears that all 
of these parameters may be treated as independent quantities in calcu- 
lating unsymmetrical landing loads. 

With regard to the wing lift at contact, both past and more recent 
studies have shown that, for civil transport airplanes, the most probable 
value of wing lift at contact is 1 g and that in 95 percent of the landings 
the wing lift was between 0.9g add 1.1 g. 

The statistical analysis indicated no significant correlations 
between wing lift and vertical velocity. 

Although a fairly substantial amount of statistical data has been 
obtained on these Initial contact conditions, the sizes of the samples 
are st il l insufficient to permit the final resolution of the effects of 
such variables as airplane type, size, wing loading, and so forth. 
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Reactions and response .- On the subject of landing-gear reactions 
and airplane response , the first problem to be discussed is the deter- 
mination of the vertical loads when the initial contact conditions are 
given. 

In recent years analytical methods have been developed for predicting 
the behavior of the landing gear during impact and the loads applied to 
the airplane. Good agreement has been obtained between the theoretical 
results and data from landing-gear drop tests. The analytical method alBO 
appears to be reasonably satisfactory when compared with flight-test data. 

Figure 5 shows a comparison of the max imum vertical loads obtained 
from flight tests of the SNJ airplane and calculated maximum loads based 
on the assumption that the airplane is a rigid body. The ordinate is the 
maximum center-of -gravity acceleration obtained for each impact, and the 
abscissa is the vertical velocity at contact. There Is considerable scat- 
ter in the flight data, about 30 percent above 4 feet per second, reflec- 
ting the effects of such factors as shock-strut binding due to drag loads, 
unsymmetrical landings, side force due to yawed landings, and so forth, 
all of which can modify the vertical loads developed. At the lower values 
of vertical velocity the experimental data are somewhat higher than the 
calculated curve, largely because of the effects of shock-strut preloading 
and breakout friction which were not considered in these calculations. 
Above a vertical velocity of about 4 feet per second, the calculated curve 
seems to lie fairly close to the mean of the experimental data. 

Figure 6 shows a similar comparison of measured and calculated maxi- 
mum vertical loads for the larger and more flexible B-29 airplane. The 
ordinate is the maximum landing-gear vert i cad-load factor and the abscissa 
Is the vertical velocity. The dashed line represents the calculated 
values. The break in the curve results from the effects of shock-strut 
breakout, which were included in these calculations . Although there is 
a scarcity of test points at the higher values of vertical velocity, it 
does appear that the calculated results are in fairly good agreement with 
the flight -test data. 

With regard to the prediction of spin-up drag loads, the most impor- 
tant unknowns have been the magnitude and variation of the coefficient of 
friction between the tire and runway during the impact. Figure 7 sum- 
marizes the available experimental Information on the variation of coef- 
ficient of friction during the wheel spin-up process. The ordinate is 
the coefficient of friction, and the abscissa is the skidding velocity, 
which Is defined as the instantaneous difference between the forward 
speed of the airplane and the peripheral speed of the tire. The shaded 
bands are envelopes of the data obtained in several different test pro- 
grams. The two relatively narrow bands represent tests with an SRJ landing 
gear in the Langley impact basin. The shorter of these two bands is for 
true forward-speed impacts; the longer one__represents impacts with forward 
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speed in combination with reverse wheel prerotation (to simulate higher 
horizontal velocities beyond the range of the impact -basin equipment). 

Hie wider dotted band shows results obtained with the B-29 airplane in 
flight tests. The lowest band represents data recently obtained with a 
B-57 landing gear at the new landing loads track. 

All these bands illustrate the fairly well-known fact that the coef- 
ficient of friction decreases with increasing skidding velocity. The 
width of the bands represents the effects of other factors which influ- 
ence the coefficient of friction, such as the slip ratio, the vertical 
load, variable-heating and contamination effects, and variations in the 
runway surface conditions, plus, of course, experimental errors. It 
will be noted that the SNJ and B-29 data, which are for relatively low 
tire pressures, are generally in good agreement with one another, whereas 
the results for the B-57 landing gear, which has a much higher tire pres- 
sure, indicate appreciably lower coefficients of friction. These results 
appear to verify a previously suspected trend, namely, that the coeffi- 
cient of friction decreases at higher tire pressures. 

Although the shapes of the bands are fairly consistent for the 
various tests, it is not as yet possible to take into account quantita- 
tively all the factors that contribute to the spread of the values of 
the coefficient of friction. Therefore, it seems that the most practi- 
cal approach at present Is to select a shape for the curve, and then 
normalize it on the basis of the maximum value. From inspection of the 
data it appears that a reasonable range of the maximums might be from 
about 0.4 to 0.9, the value to be used depending on the conditions 
involved. 


Taxiing Loads 

On the subject of taxiing loads, essentially two types of data are 
available at present; that is, measurements of runway profiles and center- 
of -gravity acceleration data from the impact portion of VGH records. 

Two possible ways of utilizing runway profile data are: (l) as 

specific displacement inputs in analog computations, or (2) as statisti- 
cal power spectra in generalized harmonic analysis. Both of these 
approaches present certain practical diff Iculties . In the case of the 
analog method, one of the major problems is to define a representative 
runway. As to the use of generalized harmonic analysis, the main stumb- 
ling blocks are as yet unsolved problems of nonlinear landing-gear 
response and the determination of a unique transfer function. 

With regard to the VGH acceleration data, several general relation- 
ships have been found which appear to be applicable to a number of 
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different types of airplanes and which, provide a common basis for esti- 
mating the loads experience in taxiing. [These are shown in the following 
table : 




AIRPLANE 

TYPES 

TAXIING TIME 

APPROX. 5 MIN 

3 

C.G. R ESPONSE 

1 TO 2 CPS 

5 

SPEED AT An MAX 

< 20 MPH 

4 

SCATTER AT CUM. 
FREQ.= I0 4 

1.5 : l 

4 


First , an analysis of the time spent in taxiing, for three transport 
airplane types, showed that this time was remarkably constant, averaging 
within a few seconds of 5 minutes. 

Second, it was found that the predominant frequency in the center- 
of -gravity accelerations - that is, the frequency at which the airplanes 
respond most to the runway roughness inputs - is also remarkably constant 
for f ive ' different types of airplanes ranging from the B -3 6 bomber to a 
small jet fighter. This predominant frequency varied only from about 1 
to 2 cps, with an average value of about 1.75 cps, despite a very wide 
range of wing bending frequencies for the airplanes involved. This fre- 
quency of about 1.75 cps appears to be associated with the natural vibra- 
tion of the airplane on its tires. 

Third, it was found that, for four transport types. in several hundred 
taxi runs, the maximum taxiing loads occurred at speeds usually below 
20 nph. 

Fourth, detailed analysis of all the loads experienced in taxiing 
for four airplane types showed that the probability distributions were 
very similar, and the scatter in the load level at a given probability 
was small, being about 1.5 to 1, or less than the scatter for such param- 
eters as vertical velocity or horizontal velocity. Consequently, it may 
be assumed that a combined distribution of taxiing loads can be used to 
represent all taxiing operations. 

Figure 8 Bhcws distributions of center-of -gravity accelerations in 
taxiing, per 1,000 flights. The ordinate is the number of times a given 
value of center- of -gravity acceleration is equaled or exceeded in 
1,000 flights; the abscissa is the center -of -gravity incremental accel- 
eration. The ordinate scale can be ratioed in direct proportion to the 
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number of flights considered. In calculating these curves, the approach 
used was first to construct a cumulative frequency distribution of all 
the individual incremental load peaks, both positive and negative, based 
on a very detailed analysis of a relatively few VGH records; these results 
were ratioed up to represent 1,000 flights, and are shown by the upper 
curve. Then, on the basis of data for 7*000 flights, a cumulative fre- 
quency distribution was obtained for the maximum taxiing load in each 
flight. This distribution of maximum loads was also ratioed to represent 
1,000 flights and is shown by the lower curve. In the limit, the two 
curves should coincide at a cumulative frequency of 1. In view of the 
nature of the assumptions involved and the fact that the two distribu- 
tions were obtained independently, the convergence appears to be quite 
good. 


If the predominant frequency of the input to the airplane is assumed 
to be 1.75 cps, one can make a straightforward response calculation for 
steady-state forced vibration and, thus, dete rmin e the variation of the 
acceleration along the span, which can then be compared with the accelera- 
tion at the center of gravity. The ratios resulting from such a compari- 
son would permit converting the abscissa scale in figure 8 from center- 
of -gravity acceleration to acceleration at any point along the span. 

As an example, such a response calculation has been made for the 
B-29 airplane; two flexible symmetrical bending modes were considered 
and a steady-state sinusoidal Input force at 1.75 cps was assumed to act 
at the landing-gear attachment point. Figure 9 shows a plot of the accel- 
eration response of the B-29 as a function of position along the span. 

Hie dashed line is the rigid-body response; the solid line is the total 
dynamic response. For the input frequency of 1.75 cps, the ratio of 
incremental wing-tip acceleration to incremental cent er-of -gravity accel- 
eration is about 7 . 9 . It is of Interest to note that data obtained in 
taxiing tests of a B-29 at the Wright Air Development Center showed 
acceleration ratios of the same order. However, there are still some 
aspects of the problem which are not fully understood; for example, even 
though the B-29 records showed frequencies of about 1.5 to 2.0 cps at 
the center of gravity, the frequencies at the wing tip were, in seme 
cases, considerably higher. 


Braking Loads 

The last topic to be discussed is the subject of braking friction. 
Recent KACA teBts with a special tow cart and with a C-12JB airplane 
have provided the data shown in figure 10. Here the maximum coefficient 
of friction in braking on dry concrete is plotted against the horizontal 
velocity. The circles are data obtained with the C-12JB airplane; the 
squares are results obtained with the tow cart at lower speeds. Although 
this type of information can be useful in the calculation of maximum 
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braking loads , no statistical data are as yet available regarding pilot 
operating practice in applying the bra k es, that is, the magnitude and 
frequency distributions of braking loads. 


CONCLUDING REMARKS 


The foregoing material summarizes the more important NACA data 
applicable to the calculation of repeated ground loads. Additional work 
will be necessary to extend the size and scope of the statistical samples 
and to fill in the several analytical gaps which still exist. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 1957* 
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VERTICAL VELOCITY FOR CIVIL AND MILITARY AIRPLANES 


PROBABILITY 



Figure 1 


AIRSPEED AT CONTACT FOR CIVIL AND MILITARY AIRPLANES 



Figure 2 
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BANK ANGLES FOR SMALL AND LARGE AIRPLANES 


PROBABILITY 



Figure 3 


ROLLING VELOCITY FOR SMALL AND LARGE AIRPLANES 


PROBABILITY 
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Figure 4 
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CALCULATED AND EXPERIMENTAL MAXIMUM VERTICAL 
LOADS IN FLIGHT TESTS OF SNJ AIRPLANE 



Figure 5 


CALCULATED AND EXPERIMENTAL MAXIMUM VERTICAL LOADS 
IN FLIGHT TESTS OF B-29 AIRPLANE 
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Figure 6 
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COEFRCIENT OF FRICTION DURING SPIN-UP 
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Figure 7 


DISTRIBUTIONS OF TAXIING LOADS PER 1,000 FLIGHTS 
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Figure 8 



NACA aw 4008 


15 


RESPONSE OF B-29 AIRPLANE 

f * 1.75 CPS 
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Figure 9 


BRAKING TESTS ON DRY CONCRETE 
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